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Q\ • Abstract 

The Bethe-Salpeter equation (BSE) for bound states in scalar theories is re- 
^ ■ formulated and solved in terms of a generalized spectral representation directly in 

Minkowski space. This differs from the conventional approach, where the BSE is 
^ . solved in Euclidean space after a Wick rotation. For all but the lowest-order (i.e., 

ladder) approximation to the scattering kernel, the naive Wick rotation is invalid. 
Our approach generates the vertex function and Bethe-Salpeter amplitude for the 
entire allowed range of momenta, whereas these cannot be directly obtained from 
the Euclidean space solution. Our method is quite general and can be applied even 
in cases where the Wick rotation is not possible. 
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1 Introduction 

The Bethe-Salpeter Equation (BSE) [|l| describes the 2-body component of bound- 
state structure relativistically and in the language of Quantum Field Theory (for an 
extensive review, see Ref. 0; Ref. is an exhaustive list of BSE literature prior to 
1988). It has applications in, for example, calculation of electromagnetic form factors of 
2-body bound states and relativistic 2-body bound state spectra and wavefunctions. 

BSEs have been solved analytically for separable kernels and for scattering kernels in 
the ladder approximation. Solutions for BSEs have also been obtained for QCD-based 
models of meson structure in Euclidean space; these solutions must be analytically con- 
tinued back to Minkowski space. It is important to note that analytical continuation back 
to Minkowski space from the Euclidean space solution is quite difficult even for the simple 
case of constituents interacting via simple particle exchange in the ladder approximation 
to the scattering kernel. In particular, BS amplitudes with timelike momenta cannot be 
unambiguously obtained from the Euclidean space solution without solving further sin- 
gular integral equations. Furthermore, for any BSE with a non-ladder scattering kernel 
and/or with dressed propagators for the constituent particles, the proper implementa- 
tion of this procedure (known as the Wick rotation |Q|) itself is highly non-trivial. For 
these two reasons the direct solution of the Minkowski space BSE is preferable. Here we 
outline such a method for scalar theories, based on the Perturbation Theoretic Integral 
Representation (PTIR) of Nakanishi |^. 

The PTIR is a generalisation of the spectral representation for 2-point Green's func- 
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tions to ri-point functions. A general n-point function may be written as an integral over 
a weight distribution, which contains contributions from graphs at all orders in pertur- 
bation theory. As any graph with n fixed external legs can be written in PTIR form, this 
must also be true of any sum of such graphs. Hence the PTIR for a particular renor- 
malised n-point function is an integral representation of the corresponding infinite sum 
of Feynman graphs with n fixed external legs. 

The scalar-scalar BSE has been solved numerically in the ladder approximation after 
the Wick rotation by Linden and Mitter Here we present Minkowski space solutions 
to the ladder BSE, which will act as a check of our implementation of the approach to 
be used here. Our numerical solutions are obtained by using the PTIR to transform 
the equation for the proper bound-state vertex, which is an integral equation involving 
complex distributions, into a real integral equation. This equation may then be solved 



numerically for an arbitrary scattering kernel ITJ]. We will restrict our consideration 
of explicit numerical solutions to the ladder approximation, although the approach is 
a completely general one. Calculations for non-ladder kernels are under way and these 
results will be presented elsewhere [|T2|. As a specific example of a scalar theory to which 
our formalism may be applied, consider the cfP'a model, which has a Lagrangian density 



where g is the (p-a coupling constant. 
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2 Formalism and PTIR 



The Bethe-Salpeter equation in momentum space for a scalar-scalar bound state with 
scalar exchange is 



where $ is the Bethe-Salpeter (BS) amplitude, and where K is the scattering kernel, 
which contains information about the interactions between the constituents of the bound 
state. 

In Eq. (^, Pi is the four-momentum of the i*^ constituent. We also define p = 
V2P1 — V1P2, which is the relative four- momentum of the two constituents, and P = pi +p2 
is the total four-momentum of the bound-state. The real positive numbers T]i are arbitrary, 
with the only constraint being that r^i +772 = 1. In the case where the scalar constituents 
have equal mass, it is convenient to choose rji = ^ = r]2, and so henceforth these values 
of the rji will be used. This notation is used, for example, by Itzykson and Zuber ||^. 

The quantities D{pi) are the propagators for the scalar constituents. We will use 
free propagators here for simplicity, although we could include arbitrary nonperturbative 
propagators by making use of their spectral representation 



where (m + /i) is the invariant mass of the first threshold m D. It is relatively straight- 
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(3) 



forward to generalise the discussion below to include p^{a) 7^ 0. We now redefine the 
scattering kernel K such that K{p,q;P) = il{p,q; P), and rewrite the momenta of the 
constituents in terms of the relative momentum p and the bound-state momentum P. 
The Bethe-Salpeter equation becomes 

^p,P) = D{^+p)D{^-p) J-^^^q,P)I{p,q;P). (4) 

I is the scattering kernel as referred to by Nakanishi in his review article 0. 

In order to convert the BSE into a real integral equation, we will need to use the 
Perturbation Theoretic Integral Representation (PTIR) for the proper bound-state vertex 



and the scattering kernel [|11|]. The bound-state vertex may be represented as 



TV'-\q,P) = yr{A-\P)q) Hda f dz- . , p., (5) 

Jo J-i [a - [q^ + zq ■ P + ^) - ze]" 

The weight function p]l^{a,z) of the vertex has support only for a finite region of the 
space spanned by the parameters a and z. A lower bound on the support is given by 
z) = unless 



a > max 



^—{m + pf + -—- (m + /i - VP2) 



[m + jjL — vP^] H [m + jj.) 



(6) 



The PTIR for the vertex (3-point) function was originally derived by Nakanishi [Q; here 
we have also assumed that the vertex for bound-states with non-zero angular momentum 
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is given by the s-wave vertex multiplied by the appropriate solid harmonic . The solid 



harmonic is an /* order polynomial of its arguments, and can be written as 



yrip) = \p\'Yrip), 



(7) 



with Y{^{p) being the ordinary spherical harmonic for angular momentum quantum num- 
bers / and m and where p = p/\p\ is a unit vector. 

We have introduced a dummy parameter n, which will be of use in our numerical work 
since larger values of n produce smoother weight functions. The fact that n is arbitrary 
can be seen by integrating by parts with respect to a; in this way weight functions for 
different values of n may be connected [Til . 

We may use the PTIR for the bound-state vertex to derive the PTIR for the Bethe- 
Salpeter amplitude, since the two are related via 



We proceed by absorbing the two free scalar propagators into the expression for the 
vertex, Eq. (^), using the Feynman parametrisation [T^. After some algebra we obtain 







+ a — (j)"^ + zp ■ P + ^) — ie 



1 n+2 ' 



(9) 
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where the weight function for the BS amphtude, <y9[^](a,z), vanishes when 



pa- 
ct < min 0, (m + /X + V P^)^ — H — — 



(10) 



To include the most general form of the scattering kernel in our derivation, we use the 



where the kernel parameters {etch; • • • , fch} are linear combinations of the integration 
variables {,^1, . . . Here we have defined, similarly to before, q = {qi — q'2)/2. The 

expression for the kernel contains a sum over three different channels, labelled by st, tu 
and us. 

The support properties of the kernel weight functions in each channel, pch, have been 
derived by Nakanishi . They will not concern us here in our pure and generalised ladder 
treatments, since in both these cases the kernel weight functions are simply products of 
delta functions. 

3 Derivation of Equations for Scalar Models 

Armed with Eqs. (j^), and (pUD , we can now derive real integral equations for 
the weight functions of the bound-state vertex and BS amplitude, both of which will be 
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PTIR for the kernel: 




(11) 



7 - («chg^ + hchP ■ q + CchP^ + 4h-P^ + Cchg ■ P + fchP ■ P) -ie 



solved numerically by iteration. 

If we consider Eq. (H), we note that it is necessary to combine four factors using 
Feynman parametrisation. Having done this, we use the so-called "self-reproducing" 
property of the solid harmonics P,p!T[]: 



/ dqF{\<^')yr{q + p) = yr{p) J dqF{\q\'). (12) 

This allows us to perform the integral over the loop momentum q. Note that in order for 
the loop-momentum integral to converge, we have the following restriction on the dummy 
parameter n: 

/<2n + 2. (13) 

This arises from a simple power-counting argument, and tells us the values of n that 
are valid for a particular partial wave. For example, for s-wave solutions we may choose 
arbitrary positive n. 

After the loop momentum integration and some algebraic manipulation, we obtain 
the result 

ip]l^{a,z) = / da dz ip^J^{a,z)]C^^{a,z;a,z). (14) 



We will omit the explicit expression for the kernel function /C for the sake of brevity. 

We can derive an equivalent equation to Eq. ([I^) , which has some advantages for 
numerical solutions of the BSE. Inserting the PTIR for the vertex, Eq. (BI), into the 
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vertex BSE 



T{p,P) = I 




-D{- + q)D{--q)T{q,P)I{p,q-P) 



(15) 



(2vr) 



and once again using the self-reproducing property of the sohd harmonics to perform the 
loop integral, we obtain the result 



can evaluate both the BS amplitude and the vertex for any momenta p and P without 
solving any additional integral equation. It should be also mentioned that our approach 
is completely independent of the choice of inertial frame, so that no Lorentz boost is 
necessary to obtain the BS vertex for a moving bound state. These are very useful 
properties when applying the BS vertex to calculations of physical processes involving 
bound states. 

4 Numerical Solutions and Results 

For our numerical study, we have specialised to the case of s-wave (/ = 0) bound 
states interacting via a pure ladder kernel. 




(16) 



Once again the explicit form of JC |]T2| will be omitted for brevity. 



Once we solve for either (/^[^'(a, z) for the BS amplitude or p]l^{a, z) for the vertex, we 



J(p,g;P) 




2 - {jp-qf 



— le 



(17) 
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which corresponds to the following fixed set of kernel parameters: a^f = 1 = Cgt, bgt = —2, 
7 = for the st-channel, with the weight functions for the tu and us channels vanishing. 
We have also solved the BSE for a sum of the pure ladder kernel and a generalised ladder 



term, which has non- vanishing fixed values for the kernel parameters dch, Cch and fch |11 • 



For numerical solution of Eqs. (0) and (|T^), it is convenient to define an "eigenvalue" 
parameter A, with A being defined in terms of the coupling as A = q'^/^Att)'^, and writing 
the kernel functions /C and K. as XK and XK, respectively. This having been done, we 
solve the BSE with the new kernel functions K and K as an eigenvalue equation, by 
iteration. Our approach enables us to obtain the eigenvalue as a function of the bound- 
state mass squared, as well as the weight functions for the BS amplitude and the vertex. 



The BSE for the amplitude has been solved in a previous work [|Tl|], with accuracy of 
the order of a few parts in 100. We have also solved here the BSE for the vertex, which 
has a simpler structure. These vertex solutions, obtained for n = 2 using an optimised 
grid and sophisticated integrator and interpolator, have an accuracy of approximately 1 
part in 10"^. Our results for the eigenvalue for the case where mcr/m = 0.5, rricr being 
the mass of the exchanged cr-particle, are shown in Table 1, and are plotted in Fig. 1. 
Some examples of vertex weight functions for various values of the bound state mass are 
shown in Fig. 2. What is actually plotted for convenience is the rescaled weight function, 
p = p]^ya'^, where here n = 2 and i = 0. The parameter t] represents the "fraction of 
binding", rj = M/2m, M = \fP^ being the bound state mass. 
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Xe 


Am 


V 


Xe 


Am 





2.5658 


2.56615 


0.9 


1.0349 


1.03497 


0.2 


2.4984 


2.49883 


0.95 




0.79528 


0.4 


2.2933 


2.29371 


0.99 


0.5167 


0.51684 


0.6 


1.9398 


1.94020 


0.999 


0.3852 


0.38530 


0.8 


1.4055 


1.40560 


1 


0.3296 





Table 1: Eigenvalues for the Wick-rotated (A^;) and Minkowski-space 
(Am) solutions of the Bethe-Salpeter Equation. The Wick-rotated values 
are from Linden and Mitter IHl. 



5 Conclusions and Outlook 

We have obtained numerical solutions of the Minkowski space BSE for scalar-scalar 
bound states in the pure ladder model. Our results agree exactly with those obtained in 
the Euclidean space treatment of Linden and Mitter 0. Our technique can be generalised 
to arbitrary scalar-scalar bound states, given that we know expressions for the kernel 
weight functions pch and the propagator spectral function p{s). The key to our approach 
is to convert the BSE from an integral equation involving complex distributions into one 
involving functions, which is numerically soluble. 

Our numerical solutions have been obtained for some simple choices of kernel weight 
function. It remains for us to carry out systematic studies for non-ladder kernels, orbital 
excitations, and comparison with other (approximate) methods such as solutions for 



separable kernels |T^. It will also be desirable to include spectral functions for the 
constituent particle propagators in a more sophisticated treatment e.g., the "dressed 
ladder" kernel, which involves the simultaneous solution of self-energy Dyson-Schwinger 
equations and the Bethe-Salpeter equation. 
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To date the spinor-spinor and scalar-spinor BSEs have been solved only in the ladder 



approximation |13|. It will be important to attempt to formulate an approach for fermions 
similar to the one outlined here, so that we may solve the BSE involving spinors for more 
"realistic" scattering kernels. We anticipate that the most challenging aspect of this will 
be generalising the PTIR to particles with non-zero spin. Since we would like to study 
the bound-state problem in QCD, for example the fermion-antifermion BSE for mesons, 
we also need to resolve the problem of incorporating confinement into the PTIR, and of 
including derivative couplings. 
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CAPTIONS FOR FIGURES 

Figure 1: The bound state spectrum for equal constituent masses and an exchange mass 
of rria = 0.5m. 

Figure 2: The rescaled weight function p = p^2^ /a^ of the bound-state vertex for 
various values of r] in the pure ladder limit, = corresponds to the case of a massless 
(i.e., Goldstone-like) bound state. 
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Eigenvalue Spectrum 




.0 0.2 0.4 0.6 0.8 1.0 1.2 
^(PV2m) 



15 






16 



